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We propose a method for the experimental generation of two different families of bound entangled 
states of three qubits. Our method is based on the explicit construction of a quantum network that 
produces a purification of the desired state. We also suggest a route for the experimental detection 
of bound entanglement, by employing a witness operator plus a test of the positivity of the partial 



transposes. 
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1. INTRODUCTION 



Entanglement, one of the central themes in quantum information processing, is well understood in low-dimensional 
systems. In dimensions 2x2 and 2 x 3, a necessary and sufficient condition for entanglement exists: the partial 
transposition test However, the properties of entanglement are much less clear in higher-dimensional systems, 

I , for which only sufficient conditions for a density matrix to be entangled are known 0, 0, • There exist higher- 
' dimensional states that, although entangled, have a positive partial transpose (PPT). Due to this property, it is 
not possible to distill any entanglement from them with local operations and classical communication (LOCC). 
Undistillable states are also called bound entangled 0. For systems consisting of more than two parties, a state may 
be undistillable even if some of the partial transposes are non-positive 7]. Even for bipartite systems, bound entangled 
states with non-positive partial transpose probably exist 8] . 

Apart from the interesting fundamental nature of bound entangled states, their usefulness for quantum information 
processing has been studied: bound entangled states can activate the distillability of one copy of a bipartite state 
with non-positive partial transpose 0, Hfl | . It has also recently been shown that one can extract a secure key from 
O |' bound entangled states fll ] . In the context of key creation, results from quantum information theory, with special 
use of bound entangled states, have recently been proven to be fruitful for insights into open classical information 
theoretical issues ^2(. Various classes of entangled states have been constructed theoretically [l^. However, the topic 
of generating bound entanglement in the laboratory and proving the produced state to be bound entangled has not 
been addressed so far. 

How does one generate a certain bound entangled state experimentally? A solution that is straightforward from 
a theoretical point of view is to consider the spectral decomposition of the state and to compose a mixed density 
matrix by creating the eigenvectors with probabilities that are specified by the according eigenvalues. However, this 
is, in general, a demanding experimental task, as one would need a source that can emit various types of product 
vectors and entangled vectors with high fidelities and well-specified probabilities. A more satisfactory approach is 
to deterministically generate a state that is the purification of the wanted bound entangled state in some higher- 
dimensional Hilbert space. The additional dimensions are provided by ancilla systems. Then, by tracing out the 
ancilla (i.e. experimentally simply ignoring the ancilla part), one arrives at the desired bound entangled state. 

In this paper we develop the latter method. Namely, we explicitly construct quantum networks that generate the 
two families of bound entangled states of three qubits introduced in [14J and pj . The first family is PPT with respect 
to any of the subsystems but nevertheless entangled, while the second family has a parameter range in which it is NPT 
only with respect to one subsystem which is not sufficient for distillation of a singlet between any two of the parties 
■ The properties of the latter states have been used recently in the context of quantum cryptography to show that 
so-called bound information exists |l2| . The networks in both cases act on a six-qubit register that is initially in state 
| 000000), and from which they generate a six-qubit pure state, such that the reduced density operator Abound of the 
first three qubits is the desired bound entangled state. 

The network for the family HI requires only eight two-qubit gates and one Toffoli gate with three control qubits, 
while the network for the family |7j requires six CNOT gates, one control-U with two control qubits and one Toffoli gate 
with three control qubits. The number of qubits and number of gates is in foreseeable reach of quantum information 
technology: at present, with NMR techniques an order-finding algorithm has been performed with 5 qubits and 6 
control gates 15]. In ion traps, 6 qubits could be provided, and control gates and simple algorithms have been 
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demonstrated [16(. 

The second step for the experimental generation of bound entangled states is to show that the generated states 
indeed carry bound entanglement. For the family of bound entangled states in [13 we discuss this issue explicitly. 
The entanglement of the state can be proved by using an entanglement witness |5(. We construct an appropriate 
witness, and provide its local decomposition which requires only four measurements settings. Furthermore, this 
family of states has a PPT with respect to any subsystem. For the experimental proof of this fact we compare three 
methods: we consider the full state estimation of the produced state pbound, the more direct spectrum estimation 
of p' hound = SPA(pb oun< i), where SPA is the LOCC version of the structural physical approximation to the partial 
transpose [rj, and finally the spectrum estimation of the partial transpose of pbound via the LOCC version of the 
network introduced in |18j . 

The paper is organized as follows. In Sec. 2 we will introduce the network that generates the class of bound 
entangled states described in 0]. In Sec. 3 we construct the entanglement witness that detects entanglement in the 
density matrix. In Sec. 4 we discuss the three different approaches to check the positivity of the partial transpositions 
of the density matrix with respect to any of the three subsystems. In Sec. 5, we construct a network that generates 
the family of bound entangled states of Ref. Q and discuss how the methods applied in Sees. 3 and 4 could be used to 
experimentally prove the existence of bound entanglement in this case. Finally, in Sec. 6 we conclude with a summary 
of our results. 



2. GENERATION OF BOUND ENTANGLED STATES 



In this section we explicitly construct the quantum network that generates the following class of bound entangled 
states [l4j| : 



21 GHZ)(GHZ \+a\ 001}(001 1+61 010) (010 |+c| Oil) (Oil |+-| 100)(100 | + -| 101)(101 1 + -| 110)(110 |Y 

c b a J 

(1) 

where | GHZ) = (| 000) + |lll))/\/2, the coefficients fulfill a,b,c > and ab ^ c, while the normalization reads 
N = 2 + a + 6 + c+ l/a + l/6+l/c. This mixed state can be generated deterministically by a quantum network 
that uses a register with three qubits plus three auxiliary qubits, all initialized at | 0) , and generates a pure states of 
6 qubits, such that the reduced density operator of the three qubits of interest is pbound- 

The procedure to generate the bound entangled state consists of two parts: a preparation stage for the first three 
qubits, and a purification stage where from the prepared state and an ancilla state a purification of Abound is generated. 
In the preparation stage one starts with the three-qubit state | 000) , and prepares a three-qubit pure state of the form 



bound/ — 



-L(| 000) + VE\ 001) + Vb\ 010) + Vc| on) + 4=1 10 °) + 4=1 101 ) + 4=1 u °) + 1 m ))- 

y/N V c Vb V a 



(2) 



This is achieved by applying certain local rotations (LU) on the three qubits, a control-U gate CU( 31 ) between qubit 
3 and qubit 1 (qubit 3 acts as the control qubit), a control-U gate CU( 3i2 ) between qubit 3 and qubit 2 (qubit 3 acts 
as the control qubit) and a CNOT gate between qubit 1 and qubit 3 (qubit 1 acts as the control qubit). This sequence 
of gates is illustrated in the left part of Fig. 
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FIG. 1: The network for creating the bound entangled state given in Eq. Q. 
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The specific form of these gates is given by 

1 1/Vb\ 

cu (3>1) = %, 2 )^ioxoi 3+ ^3( ( ( ^ + ^ ) w 

CU C3>2) = V^lOXOla + I^^^^^) ( ^^) 9 ®UX1|.. (5) 

where Ni = y/bj (1 + b), N2 = 1/ \/l + 5, = y/c/(\ + ac), and ./V4 = y/a/{a + c). The coefficients a and /3 depend 
on a, 6, c and must be chosen such that aN\N2 = PN3N4 and a 2 + (3 2 = 1. 

It is straightforward to confirm that this set of gates is constructed such that it performs the following sequence of 
transformations: 



! 000) 3 N, (l 0) + -L| 1>) JV a (| 0) + V6| 1>) (a| 0) + /3| 1)) 

|0) + -^|l))(|0) + ^|l))|0)+(^|0) + -J=|l))(|0) + y-|l;) I 



CU(3,1) -CU(3,2) 1 

CNOT (li3) 



bound/ 



In the second part of the network one first applies a sequence of three CNOT gates between the main and the 
auxiliary qubits: in this way each term of | Abound) is copied to the ancilla system. Here, the first, second and third 
qubits of the main system act as control qubits, and the first, second and third ancilla qubits act as target qubits, 
respectively: 

I VWmd} I 000) 3 ^2 Ts -L(| 000)| 000) + y/E\ 001)| 001) + Vb\ 010)| 010) + v^| 011)| 011) + -=| 100)| 100) 

x/N V c 

+ i=| 101)| 101) + 1=1 110)1 110) + |111)|111)). (6) 
Applying CNOT (4i5) and CNOT (4 ;6) then leads to 

— > -=(\ 000) I 000) + y/a\ 001) I 001) + Vb\ 010) | 010) + \/c\ 011) | 011) 
v N 

+ 1=1 100)| 111) + -1| 101)| 110) + -1| 110)1 101) + 1 111)1 100)). (7) 
V c x/c y/a 

Finally, one applies a 3-Toffoli gate, where the three system qubits are the control qubits and the first auxiliary qubit 
is the target. Its action is defined as 0] 

\a,b,c)\f)^\a,b,c)\a-b-c®f). (8) 
The resulting state of the total system is then 



Abound) 



-L ((I 000) + I 111)) I 000) + Va| 001) I 001) + Vb\ 010) I 010) + \fc\ 011) 1 011) + -f=\ 100) 1 111) 

+ -1| 101)| 110) + 1=| 110)| 101)) . (9) 

Tracing over the three auxiliary qubits, one obtains that the remaining state of the three system qubits is of the 
desired form of Eq. : 

Tr aux (|* bound )(* bound |) = ±-(2\GHZ){GHZ\ + a\ 001) (001 1 + b\ 010) (010 | + c\ 011) (011 1 



N 



+ -|100)(100| + ll01)(101| + -|110)(110|) =Pbound. (10) 
c a/ 
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The total quantum network that generates the bound entangled state pbound is shown in Fig. ^ 

Note that for the generation of this bound entangled state a more general version of the Toffoli gate can also be 
applied, namely | a, b, c) \ f) — ► exp[z6*(a, b, c)]| a, b, c)| a ■ b ■ c ffi /), because the extra phases cancel when one traces 
over the ancilla qubits after the Toffoli gate. This requires less elementary operations than the Toffoli gate ^j. The 
Toffoli gate with three controls can be decomposed into 13 two-qubit gates [2(j. We point out that in this paper we 
are mainly interested in providing a network for the generation of bound entanglement with a small number of gates, 
rather than in the optimization of this network, or the decomposition of the necessary gates into elementary single 
and two-qubit gates. The latter issue is discussed elsewhere in the literature |2lj . 



3. CONSTRUCTION AND DECOMPOSITION OF THE ENTANGLEMENT WITNESS 

In this section we will construct and locally decompose an entanglement witness W that allows to detect the 
entanglement of Abound with only four local measurements. An entanglement witness 0,0 is an operator with non- 
negative expectation value on separable states, and with negative expectation value on some entangled states. Thus, 
if we construct an appropriate witness for pbound and then measure its expectation value, the experimental result 

Tr (PFp bound ) <0 (11) 

indicates that Abound is entangled. Witnesses are observables with non-local eigenvectors which would be difficult 
to measure directly. However, witness operators can be decomposed locally |22| and thus be easily measured in an 
experiment. 

The state pbound that we want to detect has a positive partial transpose with respect to every subsystem, but there 
is no product vector | </>) in its range s.t. | (j>* x ) is in the range of Pbound (here X = A,B,C; the symbol Tx denotes 
partial transposition with respect to subsystem X, and * x denotes complex conjugation with respect to subsystem X). 
Such states are called bound entangled edge states || . Any PPT entangled state cannot be detected by decomposable 
witnesses, i.e. witnesses which are of the form W = P + Q Tx , where P and Q are positive operators. However, 
there arc methods for constructing witnesses that detect bound entangled edge states. We will follow the methods of 
Ref. 5] for the construction of our witness [1^, namely: 



W = W-eI, (12) 

where 

W = P + Q T A A +Q T B B +Q T C C 1 (13) 

where P denotes the projector onto the kernel of pb and Qx is the projector onto the kernel of p T B x ■ The parameter 
e is given by 

e = inf (e,f,g\W\e,f,g), (14) 

from which e > follows . We find that 

1 / A 1 / o , . , a 1 



w = - (| oooxooo I + 1 iii)(iii |) + y^{° 2 \ ioo) (ioo | + 1 on)(on |) + TT ^(l oio)(oio | + b 2 \ ioi)(ioi |) 
+ TT^(' ° 01 ^ 001 1 + fl2 l 110 >( 110 1) - [\ + + TTfe 5 + TT^l (' 000)(m 1 + 1 m )(° 00 1)- ( 15 ) 

Employing the Pauli operators a z = \ 0)(0 | — 1 1)(1 1, a x = | 0)(1 1 + | 1)(0 |, and a y — —i\ 0)(1 1 + i\ 1)(0 |, one can use 
the decomposition [24| 

| 000>(111 1 + I 111>(000 | = i(af 3 - * x *y* y - o y o x o y - a y a y a x ) = ±(*® 3 - \{o x + <j y )® 3 - \{a x - a v f 3 ). (16) 

For the last expression the local operators af and I (a x ± a y )/y2\ have to be measured. All the other projectors 

in Eq. (|15|l can be measured with a single af 3 measurement. Hence, the measurement of the witness requires only 4 
measurement settings. Using the methods of |24| this number can be proved to be optimal. On the other hand, if state 
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FIG. 2: Lower bound for e a=i)=1/ / c as a function of a: The left curve (1) is given by a 2 /{I +a 2 ), while the right curve (r) is the 
analytic minimum of e for 8 e — Of — 6 g . The maximal lower bound e max > 0.1069 is obtained for a = 0.3460, where the two 
curves meet. The result of the numerical minimization (n) is plotted on top of the two analytical curves, and equals the lower 
branch of them for all a. 



tomography is applied to confirm the positivity of the partial transposes (c.f. chapter IV.), then all measurements 
necessary for the witness with the first decomposition of Eq. I|16[l are already performed there. 

The last step on the construction of our witness is the computation of the value of e. We use the parametrization 
| e) = cos6> e | 0) + expire sin0 e | 1) and accordingly for | /) and | g). This leads to 



e = inf 

\ej,g) 



\{i C eCfC g ) 2 + (s e S f S g ) 2 ^) + — ^ (c 2 {s e CfC g ) % + {c e S f S g ) 2 ) + — ^ {i C ^ S f C af + ^ \& >eC f S ff ) 2 ) (17) 



+ TT^2 (( CeC / S f) 2 + a2 ( S e S f C a) 2 ) 



1 


1 - 


1- c 




T 




.2 



2 cos( 



+ 4>f + (f>g)C e CfCgS e SfSg) 



where c e j t9 = cos8 e j_ g and s e ,/,g = sin9 e j t g. In this equation the phases (j> e ,4>f ,<j> g appear only in the term 
cos((/> e + cf>f + <j>g). Therefore the phases can be chosen to be equal to zero, using the following argument: the term 

^2cos(</> e + 4>f + 4>g)c e CfC g s e s fSg^j in the above equation has to have a positive sign in order to minimize e. As the 
coefficients c e j t g and s e j^ g occur only quadratically in all other terms, all of them can be chosen to be positive. Then 
e is minimized for <j) e = <j)f — <p g = 0. We are thus left with 6 real parameters. If the parameters a, 6, c are determined 
by the experimental set-up, then the corresponding value of e can be obtained numerically by use of a multivariable 
minimization routine. 

If the parameters a, 6, c can be chosen freely, then it is advantageous to maximize e with respect to a, 6, c. Making 
the natural assumption that white noise is introduced in the preparation procedure of the state, i.e. p p = ppB + ^~g^^i 
the witness will detect entanglement in the state for p > 1 — 2e. Hence the tolerance of the witness to the presence 
of noise is enlarged by maximizing e. We searched for the maximum in the parameter range a = b = 1/c e]0,1[. 
(Remember from the definition of pbound in Eq. Q that one has to use the open interval here.) We obtain numerically 
that for a < a t h = 0.3460 the minimum is reached at e = a 2 / (1 + a 2 ), i.e. when the product state is one of the three 
possibilities | e, /, g) — |011),| 101), | 110). For a > ath the minimum of e is obtained when 6 e = Of = 9 g . These 
results are shown in Fig. |3 We find e^£=i/c — 0.1069 which is reached for ath = 0.3460. This is also the highest 
value obtained numerically when a, b, 1/c G ]0, 1[ without the restriction a = b = 1/c. For this choice of parameters 
the state mixed with white noise as described above is still detected for p > 0.786, i.e. more than 20% of white noise 
can be tolerated. 



4. TESTING THE POSITIVITY OF THE PARTIAL TRANSPOSE 



In this section we present three different methods to check the positivity of the partially transposed density operator 
Abound with respect to subsystem X = A,B,C. One possible option is to perform the full state estimation of 
Pbound [HI, and then to check whether all the eigenvalues of Pbound ^ or ^ = A,B,C are positive. This method 



6 



requires the estimation of (2 x 2 x 2) 2 — 1 — 63 independent parameters of the density operator. This can be achieved 
by performing 3 x 3 x 3 = 27 measurements on single copies of pbound, since one can write any three qubit state as 

p = - ^ Xi,j,kCTi® CTj ® CTfc, (18) 

i,j^k— 0,x,y,z 

where \,j,k = tr(p(Ji <g> <7j ® Ok), and cjq = I. The data from estimating \i t m,n with l,m,n = x, y, z can also be used 
to estimate Xo, m ,n, ^i.o.n and \i,m,o- Hence, only local measurements in the x,y,z directions have to be performed. 
One disadvantage of this option is the superfluous estimation of parameters of the density operator, since we are only 
interested in learning about the lowest eigenvalue of the partially transposed density operator. 

Another method for finding out whether pjTound > ^ ^ or X = A, B ,C is to start by applying the structural physical 
approximation (SPA) 17] to the partial transpose of Abound, and then to estimate the lowest eigenvalue of the 
resulting density operator. A structural physical approximation is a completely positive (CP) map, constructed from 
a positive, but not CP map, by adding white noise. The aim in constructing these approximations is to allow the 
physical implementation of maps which are useful in entanglement detection, but are non-physical. In this way one is 
able to bypass full state estimation when trying to detect the existence of entanglement in a given system, since one 
can estimate directly the relevant parameters, e.g the lowest eigenvalue. The construction of the SPA for a positive, 
but not completely positive map A, is as follows: 

\ d 4 A I ® I , d 4 X , , 

^ A ]^ = diATT^ + (1 -^ATT) [l8A] ^ (19) 

where d is the dimension of each of the two subsystems on which [I Cg) A] acts, and A is the absolute value of the most 

d 2 

negative eigenvalue obtained when [(I<8>I)(I<8> A)] acts on the maximally entangled state X^j=i \ / V d 2 . Each state 
| i) pertains to a d 2 -dimensional system, itself composed of two subsystems of dimension d. If one takes A to be the 
transposition map T, one finds X — 1/d. In the two-qubit case, one obtains that 

[I®T](p) = |/®/ + i[I®A](p), (20) 

which can be implemented as 

— i 2 

[I®T](p) = -Ai ® A 2 + -I ® a x a z Aia z a x , (21) 

where Ai(p) = 1/3 Yli= x y z a iP a ii an d ^-%{p) = 1/4 Sj=o x y z <J iP a i- Note that the map [I ® T] can be implemented 
using only LOCC. The extension of this construction to a system of three, rather than two, qubits is trivial. All we 
need to consider is the map 

i 2 

[I(g>I®T](p) = -I® A x ® A2 + -I®I®CT x cr 2 A 1 <T 2 o- x , (22) 
o o 

since the composition of a CP map with identity is still a CP map, and the construction of [I ® T] is independent of 
the existence of any additional systems. This map can again be implemented using only LOCC. 

Hence, in order to check the positivity of p Tx with X = A, B, C, it is enough to implement I ® I ® T on Pbound) 

and then estimate the lowest eigenvalue of p' hound = [I® I ® TUmjound)- The estimation of the lowest eigenvalue of p' 
can be achieved bypassing full state estimation, following |2^.l27l|. 

Consider a typical set-up for single qubit interferometry, conveniently expressed in terms of quantum gates and 
networks: Hadamard gate, phase-shift gate, Hadamard gate, and measurement in the computational basis {| 0), | 1)}. 
We modify the interferometer by inserting a controlled- V operation between the Hadamard gates, where V is the swap 
operator defined as V\ c/>)a \ ?P)b = \ <t>) A, V| (j>), \ ip). The control is on the qubit and V acts on the quantum state 
Q = pa® Pb- The interaction between the qubit and the environment g via the controlled-!^ leads to a modification of 
the observed interference pattern by the factor ve iQ = Tr [V^pa <g> ps)} — Tr [paPb]- The generalization of the swap 
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FIG. 3: Quantum network that estimates the non-linear functionals by LOCC. 



operation V to the shift operation (where V^\ fa)\ fa)---\ fa) = I fa)\ 0i) --- j faz-i), V| fa), i = 1, k), and the 
choice of g = p® k as the input state, allows us to estimate multi-copy observables, Tr [p k ], of an unknown state p [26|. 

Let us now extend this method to the LOCC scenario by constructing three local networks, one for Alice, one for 
Bob and one for Charlie, in such a way that the global network is similar to the network with the controlled-shift. 
Unfortunately, the global shift operation cannot be implemented using only LOCC. Thus, we will implement it 
indirectly, using the network depicted in Fig. |3| 

Alice, Bob and Charlie share a number of copies of the state pabc- They group them respectively into sets of k 
elements, and run the local interferometric network on their respective thirds of the state qabc — Pabc- ^or eacri 
run of the experiment, they record and communicate their results. 

The individual interference patterns Alice, Bob and Charlie record will depend only on their respective reduced 
density operators. Alice will observe the visibility va = Tr[/9^J, Bob will observe the visibility vb — Tr [p^] and 
Charlie will observe the visibility vp = Tr [po\. However, if they compare their individual observations, they will be 
able to extract information about the global density operator pabc, e.g. about 



Tr[f?; 



k i 
ABC\ 



= Tr 



pTbc (v{ k) 9 ® v c 



(fc) 



(23) 



This is because Alice, Bob and Charlie can estimate the probabilities Pj,-; that in the measurement Alice's interfering 
qubit is found in state | i)a, Bob's in state | j) b and Charlie's in state \ l)c, for I = 0, 1. These probabilities can 
be conveniently expressed as 



-Tr 



m 
Pabc 



+ (-i)V< fc) 



1 + (-l) ! Tr (p k A ) + (-lfTr {p K B ) + (-l)'Tr ( Pc ) + (-l) l+J Tr (p AB ) 



(24) 



+ (-l) l+i Tr (p k AC ) + (-ly+'Tr (p k BC ) + (- 



Tr(p$ 



From the latter equality it follows that 



K(Pabc) = P 
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(25) 



This is equivalent to Tr (p ABC ) = (a z ® a z <£> cr z ), where a z 



I o><o | 



1. Given that we are able to directly 



estimate Tr for any integer value of k, we can estimate the spectrum of pabc without resorting to a full state 

tomography. In our case pabc = Pbound and we need to estimate seven parameters, Trp^ ound with j = 2,3, ...,8. 
Together with Tr pbound = 1 they suffice to determine the eigenvalues of Abound Hi| . 

Even though the SPA option requires the estimation of just seven parameters of Abound, it has two potential 
experimental difficulties, the first one being the feasibility of implementing the SPA, and the second one the feasibility 
of implementing the quantum networks involving C — V gates. 

Finally, we have the option of directly estimating the non-linear functionals Tr\(p Tx ) k ] with k = 1,2,3,... and 
X = A,B, C, following |l8|. This scheme is a modification of the scheme presented in [26|, and can be also implemented 
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using only LOCC |2jfl. The main difference between the quantum network of [18|, when compared with [26J, is that 
the C - V [k) gates acting on the different subsystems do not all shift in the same direction, that is, all but one will 
shift in the direction 

V^\ ai )\a 2 )\a 3 )...\a k ) = | a k )\ ai )\ a 2 ) . . . | a fc _i), (26) 
while the remaining subsystem will shift in the opposite direction 



V- 1 W\a 1 )\a 2 )\a 3 )...\a k ) = | a 2 )\ a 3 ) ■ ■ ■ \ a k )\ ax). (27) 

In fact, the subsystem with respect to which we want to partially transpose our density operator, will be the 
subsystem shifted in the opposite direction. The advantage of this option, when compared to the SPA one, is that we 
do not have to implement any map on pbound before estimating the relevant non-linear functionals. Also, the quantum 
network used in both schemes has the same level of experimental difficulty. 



5. GENERATION OF THE DUR-CIRAC-TARRACH STATES 



Another interesting family of states of three qubits, bound entagled in a certain parameter range, was introduced in 
Q. In this section we will show how to produce it experimentally, with a method similar to the one described above. 
Using the notation from 7] , this family is given by: 



PDCT= E A 0l*0><*0 



E 

fc=01,10,ll 



(28) 



Here 1^^) = -^(l^i^O) ± \k\k 2 l)), where k\ and k 2 are the binary digits of k, and ki denotes the flipped ki. 

(Note that the state | ^q") in this notation corresponds to | GHZ) from above.) The normalization condition reads 

1. With the definitions A = An 



Aq^ + A + 2(Aqi + Aio + An) 



A > and 



s k 



__ f 1 if X k < A/2 



if A fc > A/2 

the following properties of the partial transposes hold : 

p Ts > 0, sio = & p TA > 0, 



A'Ol 



sn = ^ p Tc > 



(29) 



(30) 



A singlet state between two of the parties can be distilled iff the partial transposes with respect to the two parties 
are negative. For the following choice of the parameters 



A+- 1 - 
A ° ~ 3' 



— — 0; Aqx — An — 



(31) 



the corresponding state is inseparable with respect to the splitting A — (BC) but separable with respect to the other 
two splittings. Hence no singlet can be distilled between any of the parties and the state is bound entangled. However, 
when it is mixed with two states that are obtained by cyclic permutation of the parties it turns out that the mixture 
is inseparable with respect to any partition |28j . These properties were used recently to show that bound information 
exists and can be activated |l2j . 

Let us sketch how the states of Eq. (|28H could be prepared with our scheme. The density matrix is given by 

\ 



Pdct 





2 




















2 







Aii 




























Aoi 




























Aio 




























Aio 




























Aoi 




























Aii 





V 


2 




















2 



(32) 
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We start again with the state | 000) and produce the pure state 

I Vdct) = -^=(| 000) + 1 100)) + V^T(| 010) + 1 110)) + Oil) + I 111)) + 001) + 1 101)), (33) 



where 7 = J A^ + A . 

The state in Eq. (|33() is reached as follows: Start by a local rotation and a CNOT gate 

I 000) ^ |0)(a+|0) + a_|l»|0) CN ^ 2 ' 3) |0)(a+| 00) +a_| 11)) where LUi = I ® ( ° + a ~ I ® I. (34) 



a- — a+ 

By proper choice of the coefficients a± we can then reach | V'dct) with 3 local unitaries described below as follows 

^ I 0) ( 7 | 00) + 10) + y/2X^\ 11) + V^l 01)) ^ | Vdct). (35) 

Hence we have to choose the coefficients and the local unitaries LU2 such that 

a+ 1 00) + a_ 1 11) ^ a+\ <f>)\ r/>) + a_ | <^) | V x ) = 7l 00) + y/2>^\ 10) + v / 2A^I H) + V^l 01)) , (36) 

i.e. we have to find the Schmidt decomposition of the state on the RHS of the last equation. This state has the 
decomposition 



m=£^-> ^ c -{vki^ )- (37) 

ij x ' 

The Schmidt coefficients are the positive square roots of the eigenvalues of C T C, namely 



= ^(l±yjl- 4[( 7 2 + 2A i)(2A 10 + 2A n ) - ( 7v /2A^ + 2\A i A 10 ) 2 ]) . (38) 



Then the rotation is given by LU2 = I ® Vi ® ^2 = (| u+), | u— )) and V2 = (| | v—)). The vectors | u±) can 
be obtained from (C T C — a±T) \ u±) = and the vectors | u±) from (CC T — a±I) \ v±) — 0. The last local unitary is 

given by LU3 = H ® I ® I, where if = T j \ ^ ' H a d amar d gate. 

Now we add again three ancilla qubits in the state | 000) , and by using three CNOT gates the first three qubits are 
copied. This yields the state 

, < xot's fj_ ( . . ga , 1 1 nn x®2N , /t— /. ni m ® 2 i I 1 1 n \®2s 



IVtoCT>|000) ' (| 000)® 2 + 1 100)^) + VA i(| 010)® 2 + 1 110)® 2 ) 

v v2 



+VAw(| OH)^ 2 + I HI)® 2 ) + V^T(| 001)® 2 + 1 101)® 2 )) . (39) 



[/ = l ( A v^L I (4( ), 



Then we apply the unitary 



on qubit 4 iff the qubits 2 and 3 are in the state | 00) . A 2-controlled operation usually acts when both control qubits 
are in the state 1 1), but this can be changed by flipping the control qubits before and after the gate. This operations 
leads to the state 

i I 000) ( /VI 0) + "v/a+I 1»| 00) + -L 1 100) ( 0) - ^\ 1))| 00) (41) 
2 iinn\® 2 u. /TTTn nn\® 2 j.iiii\» 2 u. /777n nm\® 2 _i_ i im\® 2 



+ VAoi(| 010)® 2 + 1 110)® 2 ) + VAiod Oil)® 2 + I 111)® 2 ) + VAiid 001)® 2 + 1 101)® 2 ) (42) 

Then a 3-Toffoli gate flips qubit 4 iff the first three qubits are in the state 1 100). Finally two CNOT gates flip qubits 2 
and 3 iff the first qubits' state is 1 1). Tracing out the ancilla particles then yields puct- Summarizing, the procedure 
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|V>dct>|000) 



3 CNOT's, 2 CU, 3 Toffoli 



CNOTi, 2 , CNOTi,; 



' ^ (| 000) - 1 100)) I ooo) + \ -^{\ ooo) + 1 100)) 1 100) 



%\GHZ~)\000) + yjX+\ GHZ)\ 100) + v%i(l 010)® 2 + 1 101) 1 110)) 

^(1 on) 82 + 1 100)1 in)) + ooi) 82 + 1 no)| ioi)) ' 



which leads to 



Tr 4 



4,5,6| V'DCTM'^DCT 



Pdct- 



(43) 



The complete network is shown in Fig. The existence of bound entanglement for the choice of parameters in 
Eq. (|31[) can be proved by showing that the state has a PPT with respect to two subsystems, but not with the third. 
This can be proved experimentally by applying the methods of Sec. 4. 

Note that the method works for a general choice of the parameters for which the rank of the density matrix is full. 



Preparation of l% CT > 

I0> [U- 

I0> — ©— 1-@ 

io> ®-(& 



Preparation of purification of p D 



I0> — & 

I0> 

I0> 



O 0- 

4 



} 



FIG. 4: The network for creating the bound entangled state given in Eq. I28|l . Open circles for the control bits indicate that 
the corresponding gate acts non-trivially on the target if the control is 0, rather than 1 as usually (filled circles). 



6. CONCLUSIONS AND ACKNOWLEDGEMENT 



To summarize, we have presented a quantum network that generates bound entangled states of three qubits. 
Explicitly we have studied the production of the two families of bound entangled states that were introduced in 14] 
and [?J . Note that our method could be adapted in a straightforward way to the generation of other types of bound 
entangled states. As our networks consists of six qubits and several two-qubit gates, they go beyond present quantum 
information processing technology - however, it seems feasible to realize them in the not too distant future. 

We also discussed different methods of testing whether the produced states generated by the network are indeed 
bound entangled. Namely, we suggested to detect the entanglement via a suitable witness operator, and to confirm 
positivity of the partial transposes by either full state estimation, or spectrum estimation of the structural physical 
approximation of the partial transpose, or direct estimation of some non-linear functionals. 

We wish to thank A. Ekert, B. Englert, O. Giihne, D. Kaszlikowski, and M. Zukowski for discussions. We acknowl- 
edge support from Deutsche Forschungsgemeinschaft and the EU (QUPRODIS). P.H. acknowledges support from the 
ESF QIT short scientific visit grant. CM. A. is supported by the Fundacao para a Ciencia e Tecnologia (Portugal). 
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